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We present a method for phase and number squeezing in two-mode Bose systems using dissipation. 
The effectiveness of this method is demonstrated by considering cold Bose gases trapped in a double- 
well potential. The extension of our formalism to an optical lattice gives control of the phase 
boundaries of the steady-state phase diagram, and we discover a new phase characterized by a non- 
zero condensate fraction and thermal-like particle-number statistics. We also show how to perform 
amplitude squeezing in a single-mode system using dissipation. 

PACS numbers: 42.50. Dv, 37.25. +k, 03.75.Kk, 67.85.Hj 



I. INTRODUCTION 

In the usual situations, dissipation, caused by coupling 
to the environment, is considered to be a serious enemy 
to quantum-mechanical systems as it leads to a rapid 
decay of the coherence. Surprisingly, however, an appro- 
priately designed coupling between the system and the 
reservoir can drive the system into a given pure state 
[2-0|. This type of quantum-state engineering, driven 
by dissipation, recently has attracted considerable inter- 
est both theoretically [H-Q and experimentally [HI E| ■ 
A strong advantage of this approach is that the desired 
steady state is obtained without active control of the 
system. This should be contrasted with the standard 
approach in quantum-state engineering where dynamical 
control of the system is required (see, e.g., Ref. fl2|). 
Another advantage is that the target state can be ob- 
tained regardless of the initial state, making the state- 
engineering protocol insensitive to imperfections in the 
initial-state preparation. 

In this paper, we present a method for phase and 
number squeezing in two-state Bose systems using dis- 
sipation. Our work is motivated by the importance of 
squeezed states in matter- wave interferometry fl3l - fl9l ]. 
By using squeezed states, the performance of an inter- 
ferometer can be increased: Phase-squeezed states im- 
prove the accuracy of the readout of the phase differ- 
ence, and number-squeezed states make longer measure- 
ment times possible. The scheme presented here gives a 
way to perform both phase and number squeezing in cold 
atomic gases and provides an important building block 
for dissipation-driven quantum-state engineering. Addi- 
tionally, by extending our scheme to an optical lattice, we 
demonstrate that it can produce a new non-equilibrium 
phase characterized by a non-zero condensate fraction 
and thermal-like particle number statistics. 

Our method can be applied to any two-state Bose sys- 
tem, such as the one described by the two-site Bosc- 
Hubbard Hamiltonian H BK = -2JS X + US%. This 
Hamiltonian is often used to describe cold Bose gases 
trapped in a double-well potential. Here, J is the tun- 



neling matrix element, U is the on-site interaction, and 
we have introduced the SU(2) generators defined as 
S x = (a\ci2 + a\a{)/2, S y = — i(a\a,2 — a\ai)/2, and 
S z = (a\ai — a\a2)/2, where a,i(a\) annihilates (creates) 
an atom in mode i. 



II. TWO-MODE SQUEEZING 

First, we consider a two-state Bose system with a fixed 
number of particles. We assume that the system is cou- 
pled to an environment, leading to dissipative dynamics 
such that the time evolution of the density operator is 
governed by the master equation, 



dp 



i[H,p} + ^-(2cpc 



e'ep 



c) , 



(1) 



where c is the Lindblad, or jump, operator. We propose 
the following jump operator for the creation of phase- and 
number-squeezed states (hereafter called the squeezing 
jump operator): 

c ={a\ + al)(ai — a 2 ) + e(a\ — a\)(ai + a 2 ) 
=2(l + e)S z -2i(l-e)S y . (2) 

Here, e (— 1 < e < 1) is a parameter by which we can 
control the squeezing [2(|. The jump operator ([2]) can 
be realized in a system of trapped ultracold Bose gas 
immersed in a background Bose-Einstein condensate of a 
different species of bosonic atoms [U [22| . 

In order to understand the action of our squeezing 
jump operator, we first consider the ideal case where 
there is no Hamiltonian and the dynamics is driven by 
the dissipative terms alone. In the double- well setting, 
for example, this can be achieved by making the potential 
barrier between the wells high enough so that J = 0. The 
interaction U can be reduced using Feshbach resonances. 
We start by calculating the amount of squeezing in the 
steady state. It can be characterized by the normalized 
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FIG. 1: (Color online) Time evolution of the squeezing and 
purity for N = 100, J = and (a) e = -0.50, (b) e = 0.74. 
The solid (dashed) line shows the numerically obtained exact 
solution for [7/7 = (U/f = 0.5), and the dashed-dotted line 
gives the analytical value of £p,n evaluated using Eq. (J5J. The 
initial state is the coherent state aq |e=o- 
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respectively, where (A5 2 J = (S$ jZ ) - {S VtZ ) 2 |2J. In 
this work, the average spin is always parallel to the x 
axis, (S) = (S x )e x . Using the coherent-state approxi- 
mation (S 2 ) ~ (S^) 2 with (S x ) ~ N/2 + O(7V ), which 
holds for any value of e provided the number of par- 
ticles N 1 [HI and the truncation scheme based 
on the Bogoliubov backreaction formalism (SiSjSk) — 
(SiS^iSk) + {Si) (SjSk) + (SA) {Sj)-2(Si} <&) (5 fc > IT 
[26j . the equations of motion for (Sy) and (S z ) become 



-4/V 7 (l - e 2 )(S 2 J+/V 2 7 (l±e) 2 . (4) 



In the second term on the right-hand side, the upper 
sign corresponds to y, and the lower one corresponds to 
z. The time constant r of (Sy) and (S z ), thus, is r _1 ~ 

4JV"7(1 - e 2 ). Since (S y ) = (S z ) = throughout the time 
evolution, in the steady state, we get 

/V lie 



( AS y,z) ry,z, - 4 1=| _ eI 



(5) 



where, again, the upper sign corresponds to y and the 
lower sign corresponds to z. We see that a phase- 
squeezed state characterized by £p < 1 is obtained for 
e < 0, while e > yields a number-squeezed state for 
which £n < 1. 

In general, the steady state of the time evolution is a 
mixed state. The amount of mixedness can be quanti- 
fied using purity P[p] = {(N + l)Tr[p 2 ] - 1}/N. For a 



pure state, P[p] — 1, while the completely mixed state 
gives P[p] = 0. The solid lines in Fig. [T]show the time 
evolution of the squeezing and purity for TV = 100 (and 
J = U = 0) obtained by numerically solving the master 
equation (JXJ) with the squeezing jump operator @. We 
have set e = —0.5 (upper panels) and 0.74 (lower pan- 
els). The value e = 0.74 is chosen to demonstrate that our 
method can produce the same (and even larger) amount 
of squeezing than a different method used in an exper- 
iment with cold atomic gases [It}- This figure clearly 
demonstrates that, using the jump operator ([2]), we can 
obtain almost pure phase- and number- squeezed states. 
Note that £p and £n of the steady states are predicted 
very accurately by Eq. (0). 

We now discuss the exact form of the steady state and 
show that our jump operator |(5J) drives the system into 
a squeezed state. The jump operator can be written as 

c = A^/ee xSx S z e~ xS:c , where \ = arctanh[(l — e)/(l + e)]. 
If N is even, one of the eigenvalues is equal to zero, the 



corresponding normalized eigenstate is 



'sq 



\N/2). 



Here, we use the notation \ j) for a state that has j parti- 
cles in mode 1. State ^> sq is a stationary state of the dy- 



namics and can be written as </> sq = J2n^-N/2 a n\N/2 



n), where 
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and A is a normalization factor. State (j) sq is a phase- 
squeezed (e < 0), a number-squeezed (e > 0), or a coher- 
ent (e = 0) state. When N is large, we can approximate 



(2^(A^ 2 )) 
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(7) 



In deriving these results, N was assumed to be even. If TV 
is odd, the steady state can be a mixed state. However, 
the non-zero elements of the vector c</> sq scale as e^ JV+1 ^ 2 , 
and 4> sq becomes an approximate dark state of the jump 
operator for large N. Thus, in the large N limit, 



is 



a steady state, regardless of whether N is even or odd. 

Now, we introduce the two-site Bosc-Hubbard Hamil- 
tonian in addition to the dissipative terms. For simplic- 
ity, we assume J = 0. We have checked numerically that 
the effect of non-zero J is much less important than that 
of non-zero U and does not change the results qualita- 
tively. 

The interaction term has the Fock states \j) as the 
steady states, while the dissipative part drives the sys- 
tem toward state sq , which is close to the coherent state 
0sq|e=o- The steady state resulting from the competition 
between these terms is a mixed state. The N dependence 
of the purity and squeezing of the steady state decreases 
with increasing N and, at large TV, these quantities are 
determined only by e and J7/7- The disappearance of the 



3 



1.0 
0.5 
0.0 
-0.5 



-1.0 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

U/y 



FIG. 2: (Color online) Purity of the steady state as a function 
of e and C//7 for N = 100. 



./V dependence can be seen by noting that both the in- 
teraction energy and the dissipative part scale as ~ TV 2 . 
The dashed lines in Fig. Q] show the time evolution of the 
squeezing and purity for the non-zero interaction. We 
have chosen N = 100. At this particle number, the iV 
dependence of the purity of the steady state is already 
small for the values of e used in this figure. We see that it 
is possible to obtain phase- and number-squeezed states 
with high purity even for non-zero U, provided C//7 is 
small enough. To quantify the effectiveness of the squeez- 
ing jump operator ([2]), we plot the purity of the steady 
state at N = 100 as a function of e and C//7 in Fig. [5] 
For fixed U/j and |e|, a positive e leads to larger pu- 
rity than a negative one and allows us to obtain a high 
amount of number squeezing without losing purity even 
when U ~ 7. This is because, for positive e, the dissipa- 
tive term favors number-squeezed states that are closer to 
Fock states than phase-squeezed states. Fock states are 
eigenstates of the interaction term and, therefore, steady 
states in the absence of dissipation. 



III. SINGLE-MODE SQUEEZING 

Next, we consider single-mode squeezing using dissipa- 
tion. In the absence of a Hamiltonian, amplitude squeez- 
ing can be achieved using the jump operator, 



c = a + ea) + iX = c + iX. 



(8) 



This is equivalent to having the jump operator c and the 
Hamiltonian H = VL(c + c')/2 with f2 = 7A. As in the 
two-mode case, e G (—1,1) controls the squeezing. At 
e = 0, we obtain a coherent state with the amplitude 
(a) = — iX [27| . The equation of motion for (a) is 



d(a) 



= -J[(l-e 2 )(a)+zA(l + e)] 



(9) 



We see that the time constant of (a) is r 1 = 7(1 — e 2 )/2, 
and 



(a) 



A 



l-e ' 



(10) 



in the steady state. The equations of motion for normal- 
ordered correlation functions (a^ rn a n ) do not contain 
higher-order correlation functions and form a closed set of 
equations. Thus, the average particle numbers (a^a), the 
number fluctuation (An 2 ) = ((a^a) 2 ) — (a' a), and the 
second-order coherence function = (a^ 2 a 2 ) / (a^ a) 2 
for the steady state can be calculated exactly, 
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(11) 

(12) 
(13) 



For > e (> e*), where e* < is a value such that 



(At 



X 2 , we obtain a number-squeezed state (At 



< 



A , which has a non-classical nature, characterized by 
< 1. If e > 0, we obtain a number-anti-squeezed 



fJ 



(2) 



state (An 2 ) > A 2 . 



IV. LATTICE SYSTEM 

Finally, we consider the application of our squeezing 
jump operator ((2) to cold Bose gases in an optical lattice 
described by the Bose-Hubbard Hamiltonian. A natural 
extension of c for the lattice system is the following jump 
operator acting on sites i and j: 

Cij = 0] + - %) + - a])(a,i + aj). (14) 

Here, a; is the annihilation operator at site i. The time 
evolution is given by the master equation, 



(7/2) £> 



■ijP C i 



PCijCij) 



(15) 

To study the qualitative effect of the squeezing jump 
operator (|14[) . we employ the generalized mean- field 
Gutzwiller approach It consists of a product ansatz 

for the density operator p = <S> i f>i, where pi = Tr^j[/3] 
are the reduced local density operators for site i and the 
site-decoupling approximation H = hi, with 

hi = - J ((oi')a] + {a\')h) ~ P^i + {U/2)hi(hi - 1) 

(16) 

and hi = d|a,. 

We are interested in the region of the higher filling 
factor h = (h) > 3 where the filling factor dependence 
becomes small and we can obtain a universal result. Be- 
sides, it has been shown that distinct commensurability 
effects are absent even at low h In the following cal- 

culations, we consider a homogeneous system (hereafter, 
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FIG. 3: (Color online) Steady states of the dissipative dynam- 
ics in an optical lattice driven by the squeezing jump operator 
f| 1 4p . We have set n = 4. Panels (a) and (b) show the num- 
ber fluctuation (An 2 ) 1 / 2 and the condensate fraction \t/j\ 2 /n, 
respectively, for J/7 = 1. In (b), the blue plus sign, red 
box, and green cross correspond to the parameter values used 
in Figs. |4jb), ( c )> an d (d), respectively. Panel (c) gives the 
non-equilibrium phase diagram. The dashed lines show the 
boundary between the phases with \tp\ 2 = (thermal state) 
and \ip\ 2 0. The phase with a thermal-like condensed state 
obtained for e = —0.2 is shown by the filled yellow region. 



we omit the site index i) with n = 4 as an example. 
Because of the large fi, we expect that the generalized 
mean-field Gutzwiller approach is quantitatively reliable 
for higher dimensions. We choose a pure coherent state 
as the initial state of the local density operator and study 
the resulting steady state of the time evolution. 

Figure Eta) shows the local number fluctuation 
(An 2 ) 1//2 = ((n 2 ) — n 2 ) 1 / 2 as a function of U/zj (z is 
the coordination number) for J/7 = 1 and e = —0.2, 0, 
and 0.2. As in the two-mode case, positive (negative) e 
yields smaller (larger) (An 2 ) 1 / 2 corresponding to num- 
ber squeezing (anti-squeezing) . In Fig. E(b), we show 
the condensate fraction |i/>| 2 /n of the steady state. Here, 
tp = (a) is the order parameter. Similar to the two-mode 
system, the interaction term favors states with a definite 
number of particles and, thus, suppresses the off-diagonal 
order. Hence, \ip\ 2 /n decreases monotonically with in- 
creasing U and finally vanishes. The boundary between 
the phase with zero and non-zero \ijj\ 2 /n is shown by the 
dashed lines in Fig. EJc) for various values of J. We find 
that the region of the condensed phase \ip\ 2 /n > grows 
(shrinks) for positive (negative) e [28|. This allows us to 
control the phase boundary of the non-equilibrium phase 
diagram Fig. EK C ) by changing the value of e. 




FIG. 4: (Color online) Matrix element p m „ of the local density 
operator p of the steady state driven by the squeezing jump 
operator ()14|) with negative e. Here, we set e = —0.2, n — 4, 
and J/7 = 1. Panel (a) shows the particle-number statistics. 
Panels (b), (c), and (d) show \p m n\ of the condensed state, 
thermal-condensed state, and thermal state, respectively. 



We note that the behavior of \4>\ 2 /fi depends strongly 
on whether e < or e > [see Fig. E^b)]. In the lat- 
ter case, \ip\ 2 /n decreases suddenly to zero (solid and 
dashed lines) while in the case of negative e, \ip\ 2 /n first 
drops to a non-zero value and then gradually decreases to 
zero with further increasing U/zj (dashed-dotted line). 
To better understand the properties of the state with 
a small but non-zero \il>\ 2 /n, we show the matrix ele- 
ments p mn (in the Fock-state basis) of the steady-state 
density operator p at U/zj = 3.6 for e = —0.2 and 
J/7 = 1 in Fig. |4|c). For comparison, we also show p mn 
for the condensed state with large \4>\ 2 /n at C//27 = 3 
and a state with \i[>\ 2 /n = at U/zj = 5. The state 
with |-0| 2 /n = is well described by the thermal state 
Pnn = n n /(n + l) n+1 , which does not have off-diagonal 
elements [see Figs. SJa) and0|d)]. Figure HJa) clearly 
shows that the diagonal elements p nn at U jz'y = 3.6 are 
very close to p nn of the thermal state. Simultaneously, 
the steady state at U/z^ = 3.6 has non-zero off-diagonal 
elements [Fig. Hfc)] . We call such a state the thermal- 
condensed state. It is characterized by almost ther- 
mal particle-number statistics and a non-zero condensate 
fraction. We also note that the thermal-condensed state 
has slightly larger (An 2 ) 1 / 2 than the thermal state [see 
Fig. E(a)]. The thermal-condensed state appears as a 
new phase in the non-equilibrium steady-state phase di- 
agram for e < 0. It is represented by the yellow region 
in Fig. Etc) . This phase is separated from the ordinary 
condensed phase by a jump in the order parameter [see 
dashed-dotted line in Fig. E£b)] and it differs from the 
thermal phase by having a non-zero order parameter. 
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FIG. 5: (Color online) The ratio of the exact values of (S x ) 

and y to their approximate value, N/2. The parameters 
are the same as for the solid lines in Fig. [T] N = 100, J = 0, 
and U/j — 0. The upper panels correspond to e = —0.50, 
and the bottom ones correspond to e = 0.74. 
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FIG. 6: (Color online) The ratio of the exact values of 
(S X S%) + (S%S X ) and (S X S%) + (S 2 Z S X ) to their approximate 
values given by Eqs. (|A3() and <)A4|I . respectively. The param- 
eters are the same as for the solid lines in Fig. [1] N = 100, 
J = 0, and U/j = 0. The upper panels correspond to 
e — —0.50, and the lower ones correspond to e = 0.74. 



V. CONCLUSION 

We have proposed a way to produce number- and 
phase-squeezed states in a two-mode Bose system us- 
ing dissipation. When applied to an optical lattice, our 
scheme can be used to control the phase boundaries of 
the steady-state phase diagram. It also allows us to real- 
ize a new phase characterized by a non-zero condensate 
fraction and thermal-like particle-number statistics. 
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Appendix A: Expectation values of S x and S x and 
Bogoliubov backreaction formalism 

In the derivation of Eq. ([5]), we have used the following 
approximation: 

(S x ),^f<J%)~^. (Al) 



In order to illustrate the validity of this approximation, 
we plot the ratio of the exact and approximate expec- 
tation values in Fig. [5] The values of the parameters 
are the same as in Fig. [TJ Note that the approximation 
(|A1[) holds very well, and the relative error is less than 
1% for e = —0.5 (upper panels) and less than 2.5% for 
e = 0.74 (lower panels) throughout the time evolution. 
We have observed that the error of the approximation 
(jAip for the final steady states scales roughly as ~ 1/7V 
and, thus, becomes even smaller for larger numbers of 
particles. 

Another approximation used in our article was based 
on the so-called Bogoliubov backreaction formalism [25| , 
which approximates the expectation value of a product 
of three operators as 

-2{S i )(S j )(S k ). (A2) 

In the derivation of (A5^) in Eq. ((SJ), we used 

(SJD + (S 2 y S x ) ~2((Sj y ) + (S y S x ))(S y ) + 2(S x )(S 2 y ) 
-±{S x )(S y ) 2 , (A3) 

and, in the derivation of (AS 2 ), we used 

(S X S 2 Z ) + (S 2 S X ) c±2((Sj z ) + (Sj x ))(S z ) + 2(S X )(S 2 ) 
-A(S X }(S Z } 2 . (A4) 
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In Fig. |H1 we plot the exact values of (S x Sy) + (SyS x ) 

and (S x Sl) + (S^S X ) divided by their approximate val- 
ues given by Eqs. (|A3|) and (|A4|) . respectively. The values 
of the parameters are the same as those in Fig. [TJ We 
see that the relative error is less than 3% for e = —0.50 
(upper panels) and is less than 7% for e = 0.74 (lower 

panels). As for (S x ) and J (S%), the error of the ap- 
proximations (|A3[) and (|A4j) for the final steady states 
scales roughly as ~ l/N, which justifies the use of the 
approximation for large N as in the case of Fig. [TJ 



Finally, we note that the approximations are valid for 
the final states of the time evolution regardless of the 
initial state. In the calculations of Fig. [TJ we chose the 
coherent state as the initial state. The system is, how- 
ever, driven to the squeezed states no matter how the 
initial state is chosen [the proof has been given in the 
paragraph of Eq. ([6]) in Sec. [ITJ. The squeezed states sat- 
isfy the approximate equations (|A1[) and (|A3[) - (|A4[) to 
the accuracy shown in Figs. [5] and [51 
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